Adiabatic radiation reaction to the orbits in Kerr Spacetime 
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Geodesic motion of a point particle in Kerr geometry has three constants of motion, energy 
E, azimuthal angular momentum L, and Carter constant Q. Under the adiabatic approximation, 
radiation reaction effect is characterized by the time evolution of these constants. In this letter we 
show that the scheme to evaluate them can be dramatically simplified. 
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I/-) ■ I. INTRODUCTION 

o : 

A supermassive black hole (SMBH) accompanied by a compact object (CO) is among the most promising candidates 
for gravitational wave sources. This system may provide us the best opportunity for testing general relativity in the 
strong gravity regime. For this purpose, however, we need an accurate theoretical prediction of its waveform. 

To investigate gravitational waves from SMBH-CO binary system, we use the black hole perturbation method. The 
background geometry is Kerr spacetime, and CO is described by a point particle. In the lowest order in mass ratio, the 
particle moves along the background geodesic. In the next order its orbit deviates from the geodesic due to radiation 
' ■ reaction effects. 

In the Schwarzschild background, the particle's orbit can be characterized solely by energy E and azimuthal an- 
■ gular momentum L. We can evaluate the orbital evolution from the change rates of energy and azimuthal angular 
momentum, dE/dt and dL/dt. In the adiabatic approximation the energy and the angular momentum that a particle 
0^ . loses are equated with the ones radiated to the infinity or into the black hole horizon as gravitational waves since 
there are conservation lows for E and L including the gravitational field. In this manner we can determine dE/dt and 
dL/dt from the asymptotic behavior of the radiated gravitational waves in the Schwarzschild case. 

On the other hand, Carter constant Q is also necessary to specify a geodesic in Kerr background. Since we do 
not have a conservation law corresponding to Carter constant, the change rate of Q is not directly related to the 
asymptotic gravitational waves. Instead, we need to directly calculate the self-force acting on the particle p). Though 
the prescription to calculate the self- force is formally established LiLf. performing explicit calculation is not so straight 
forward. 

Gal'tsov employed the radiative part of the metric perturbation, which was introduced earlier by Dirac p|, to 
calculate dE/dt and dL/dt. He showed that this scheme correctly reproduces the standard results obtained by using 
the conservation laws. Recently, Mino gave a justification for applying the same scheme to dQ/dt for bound orbits [||. 
(See also Ref. Q-) However, actual implementation of dQ/dt calculation again is not so straight forward. In this 
letter we derive a rather simple new formula for the adiabatic evolution of Carter constant. 

II. BACKGROUND 

We consider the background Kerr spacetime in the Boyer-Lindquist coordinates: ds 2 = —(1 — 2 Air /E)d£ 2 — 
(4Mar sin 2 0/T)dtdxp + (E/A)dr 2 + Ed# 2 + (r 2 + a 2 + 2Ma 2 r sin 2 0/E) sin 2 Odtp 2 , where E = r 2 + a 2 cos 2 and A = 
r 2 — 2 Air + a 2 . Here Al and aAl are the mass and the angular momentum of the black hole, respectively. There are 
two Killing vectors ^ := (<9t) M and £^ := (<9 ¥ ,) /i . In addition, Kerr spacetime possesses the Killing tensor, := 

2£Z( M n„) +r 2 g fiv , where the parentheses denote symmetrization on the indices enclosed, and l 11 := (r 2 + a 2 , A, 0, a) /A 
and n M := (r 2 + a 2 , —A, 0, a) /2E are two radial null vectors. Killing tensor satisfies the equation K^.^ = 0, 

We consider motion of a point particle, z a (r) — (t z (r), r z (r), 6 z (t), 4> z {t)). Here r is the proper time along the 
orbit. For geodesic motion, there are three constants of motion, 
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Q := K ap u a uP = {L aE f n2 &z)2 + a 2 cos 2 9 Z + E V) 2 , (2.3) 

sin # z 

where u Q := dz a /dr. In addition, we define another notation for the Carter constant, C := Q — [aE — L) 2 . For orbits 
on the equatorial plane C vanishes. 



III. GEODESIC MOTION IN KERR SPACETIME 

Introducing a new parameter A defined by dX = dr/Y,, the geodesic equations become 

(£)'-*<'■>■ =«<-»».>. m 

^_^*,.,.- i) + i + £p( r .), fc = - aE+ -i_ + iip W , (3.2) 

where P(r) = £(r 2 + a 2 ) - oL, R(r) = [P(r)f - A[r 2 + Q] and 6(cos6») = C - (C + a 2 (l - E 2 ) + L 2 ) cos 2 9 + a 2 (l - 
E 2 ) cos 4 9. It should be noted that the equation for the r-component and the one for the ^-component are decoubled 
by using A. Both R(r z ) and 0(cos£? z ) are quartic functions of their arguments. Hence both solutions are given by 
elliptic functions. For bound orbits, we can systematically expand r z and cos# z in Fourier series. 
The other two equations l|3.2|) are integrated as 

t z (X) = tM (A ) + t«(A) + ^A, 

^(A) = ^')(A) + ^)(A) + ^^A, (3.3) 

where (• ■ •) means time average along the geodesic. t^(X) := J dX{(r 2 + a 2 )P{r z )/ A —{{r 2 + a 2 )P(r z )/ A)} and 
fw(A) ■= — J d\{a(aE sin 2 9 Z — L) — (a(aEsm 2 9 Z — L))}, are periodic functions with periods 2ir£l r and 2irQg, respec- 
tively. Functions ip^ and tp^ are also defined in a similar way. 



IV. ADIABATIC EVOLUTION OF CONSTANTS OF MOTION 

In Ref. it was shown that the adiabatic radiation reaction to the constants of motion P = {E,L,Q} can be 
evaluated by 

where h 1 ^ is the radiative part of the metric perturbation defined by half retarded field minus half advanced field, 
i.e., k™" :— (h 1 ^ — h^ v )/2. Radiative field is a solution of source free vacuum Einstein equation. The singular parts 
contained in both retarded and advanced fields cancel out. Therefore we can avoid the tedious issue of regularizing 
self-force. f a is a differential operator, 



It was shown by Gal'tsov Q that 



Calculation of E and L 
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where ip(x) = Y,u tl u u h llu /2 and (ut,u r ,ug,u v ) := (-E, ±y/R(r)/A, ±^Q(cos9)/ sia6,L). This vector field it M is an 
extension of the four velocity of a particle in the sense that it satisfies u M (z(A)) = u fl (X). Since in fact u r and ug, 
respectively, depend only on r and 9, we can verify the relation, Ua-p = u/3 :a - 
Furthermore, Gal'tsov has shown pJj that ip(x) is given by 

= 1 / £^ E Ci» / (4-3) 

where <^ />m (x) = Eu^u^x)^ m (x) = Eu"(r, 0)u"(r, 0)r^(r, 0) e -^*+™¥>A 2 _ 2 iC /m (r) _ 2 S wA)B (fl). 7T w/ , m is 
an appropriately normalized mode function of metric perturbations, which is constructed by applying a second rank 
differential operator t*„ to a mode function of Teukolsky equation Jg ■ The method to solve Teukolsky equation is 
well established . We can calculate the contribution from waves absorbed into a black hole by replacing the in-field 
to the up- field in Eq. (|4.3|l . 
We can express 4>w,l,m as 

Kt,m{^)) = e- l{ ^ W - m ^ W) ^r 7 XX),dr 7 JdX(X),e z (X),de z /dX(X)). (4.4) 

In the folloing text, we abbreviate dr z /dX and d9 z /dX from the arguments for brevity. Here the exponent con- 
tains t (r) (A), t^(X), ip {r) (X) and <p {e) (X). Since r z , and tp^ (9 zi and tp^) are periodic functions with 
period 271-fi" 1 (2 7 ril~ 1 ), we can expand e - 4 ^(* (r) ( A '-)+ t<< ' ) ( A e))- m ^ ( '' > ( A -)+v )<< ' ) ( A e)))$ Am (r z (A J ,), 9 z (X e )) into Fourier 
series as {dt z /dX)J2n r ,n e Z"^e m ^ A - +me0 <> Ae . Therefore we obtain exp[-iu(&)(\) + t^ 9 \X)) + im{^ r \X) + 
tpW(X))}<P t „ m (r z (X),9 z (X)) = (dt z /dX) £L 2«r.«j. e i(nrn r +n.n.)A i Using this expansion, we obtain 



/ 



dAV22m(*M) = £ 2^(c-<-^)Z^, (4.5) 



»,..?!„ 



where = (dt z /dX) 1 (m (dip z /dX) + n r Q r + n e Q e ) and Z^"" := Z^,^, Substituting Eq. (JOJ into 

Eq. 14.211 with Eq. I]4.3[l . and integrating it over u>, we obtain 



e i((n r -n' r )n T + (n e -n' g )n e )\ g-j 



n r ,riQ : n' ,71'f. 



Now integration over A is straight forward. We finally end up with the well known formula except for the overall 
normalization depending on the definition of the mode function: 

f ) = - £ wi a - (4.7) 

£,m,n r ,ng 

In a similar manner, a formula for the angular momentum loss rate is obtained as 

B. calculation of Q 

The expression for the radiation reaction to the Carter constant can be cast into a form analogous to the energy 
loss rate and the angular momentum loss rate. Substituting 

U = \{d v h a p)u a uP - ^-(h^u") - \u v ^(u a uPh a p), (4.9) 

I (XT A (IT 



the evolution of Carter constant is given by 

»/. . . 1 

(4.10) 
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Here means that terms which become a total derivative or 0(h 2 ) are neglected. Furthermore we can show that 
the second term in the above equaiton vanishes by using the facts u a; /3 = up- a and Kt^.p) = 0. Thus we obtain 



dQ 
dX 



lim — 



-T 



(4.11) 



which is written in terms of ip(x). Hence, we find that the change rate of Carter constant is obtained by replacing 
— £u\d<x m the above expression for (dE/dt) given in Eq. (|4.3|) with 2T,K^u )J, d u T l ~ 1 . Then, what we have to evaluate 
is 
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x— *z(\) 



dX 
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(4.12) 



In the last step the last term have been integrated by parts using u^d^ — E 1 d/dX, which is valid after substitution 
ofz(A). 



V. FURTHER REDUCTION 



Further simplification is possible. For an arbitrary function of r z and 9 Z , we have 



r lm- / dXe X p[-iu2> n n z (X)+im<p z (X)]f(r z (X),6 z (X)) 



x exp 



T 

dX r 



dXg exp 



-in r n r X r - iuj%' nB t( r) (\ r ) + irruf {r) {X r ) 



-in e n e X e -iuj2' ns t (e) (Xe) + imtpW (X e ) f(r z (X r ),9 z (X e )). 



(5.1) 



This relation can be easily verified by substituting exp[-iu%' n ° (t^(X r ) + t^(X e )) + im(<p^{X r ) + 
^ [e) {X B )))]f{r z {X r ),e z {X 6 )) = En r ,n e /nr,n e e^^+ne^X _ Thenj by using the above re l a tion, the A-integral in 
Eq. (|4.12|) can be decomposed. The part of A r -integral takes the following form, and it can be integrated by parts as 



dr z 
dX 

dX r 



d r exp [—in r Q r X r — iiot + imip] f(r, ( 



r-»r,(A r ),t-+t M (>»•)>¥>-► ¥> M (V) 



d fdt^X n fdJ^ 



x exp [— in r Cl r X r — iut + imip] f(r, 9) 



(5.2) 



r _ >r:s (A r ),t->.t(*-)(A».),p->»jM(A,.) 



The first term in the last line is a total derivative, and therefore does not contribute to the average over a long period 
of time. Combining the above results, we obtain 



dQ 
dt 



= 2 



dX I t-,oo T J_ T 
(r 2 + a 2 )P\ I dE 



(r 2 + a 2 )P\ /aP\ . 



tp(x) dX 

%—*z(X) 



A 



i ill \ / aP\ I dL \ \ — > 



n r ,ng m 
z,m,n r ,ne 



(5.3) 



This expression is as easy to evaluate as (dE/dt) and (dL/dt). To evaluate the last term, we have only to replace m 
in the expression for (dL/dt) with n r Q r . 
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VI. CONSISTENCY CHECK 

We know that a circular orbit stays circular under radiation reaction ^(j- This condition becomes (dQ/dt) = 
(2(r 2 + a 2 )P/A)(dE/dt) - (2aP/A) (dL/dt). For circular orbits Z™ r m n< * = when n r ^ 0. Therefore the last term in 
Eq. I|5.3[l vanishes. Then Eq. (|5.3() agrees with the above condtion that a circular orbit stays circular. 

We also know that an orbit on the equatorial plane does not leave the equatorial plane by symmetry. This can 
be clearly seen by rewriting the above formula in terms of C. An identity (dt z /d\)(dE/dt) — (d<p z /dX) (dL/dt) + 
J2e m n r n e {( n '•^'• + n s^e)/ UJ m } I -^"m™" 2 = follows from the definition of u%[> ne with the aid of the expressions 
for (dE/dt) and (dL/dt) given in Eqs. I|4.7|l and JO}. Using this identity, we have 

(f) - -><«W..>(§ )_ 2(ico ^)(§)- 2E ^ W (6 .d 

From this equation it is manifest that (dC/dt) = when 8 = tt/2. Notice that Z™^ 1 " ^ only for ne — in the case 
of equatorial orbits. 
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